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CHAPTER I 
DmODUCTION 

1.1 Motivation and Relevance of Tlieais 

llie LACATE (Lower Atmosphere Oonposition and Tenfjerature Experi- 
ment) mi,ssion was a high altitude balloon platform test which enployed 

A 

an infrared radicrveter to Si^nse vertical profiles of the concentrations 
of selected atmosEheric trace constituents and temperaturfis. Ihe con- 
stituents were measured by Inverting infrared radiance profiles of the 
earth’s horizon. Ihe radiometer line-of-sight was scanned vertically 
across the horizon at aFproxlmately 0,25°/seoond. Hie relative vertical 
positions of the data points making up tlie profile had to be determined 
to approximately 20 arc seconds. 

The balloon systan for accomplishing the mission is shewn in 
Pig. 1.1-1, It consisted of: (a) a 39 million cubic feet (zero pres- 

sure) balloon, (b) a load bar containing the balloon control equipment, 
(c) a package containing additional balloon ccntrol electronics with 
gondola reoo'/ery parachute, and (d) a gondola containing the research 
equipment. The bcdloon was designed to lift the payload to a float 
altitude of approximately 150,000 feet. 

Instrumentation to determine the attitude pf the balloon platform 
consisted of a magnetaneter and 3 orthogonally oriented precisicn rate 
gyros. Ihe three rate gyros were employed to obtain an accurate time 
history of the angular velocity carponents of the research platform for 
subsequent data reduction and attitude detemlnation. 

The main problem in the LACATE experiment is to determine the in- 
stantaneous orientation (i.e., the attitude) of the instrumentation 
platform with respect to a local vertical. Moreover, this orient(:ition 
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most be determit':^ to an accuracy of 1^. Ckice this is known # the orien- 
tation of the line-of-sight of the r'jdicp«*t#»r c«m be determined since 
its relative motion with respect to the platform 4s prescribed, 

1.2 Present Status of the Problem 

Stablizing the balloon research platform or predicting its orien- 
tation is a major problem which must be solii!^ in all balloon borne 
experiments requiring line-of-sight instrumentation. Feedback oontacol 
systems have been used to stabilize the balloon platform with respect 
to* an inertial reference frame . Stability is cAjtalned by suspending 
the platfom at its center of gravity and employing some control system, 
control system instnmentation includes sensors (rate gyros, digital 
star trackers, etc.) m'a reactiori wheels for torquing the platform. 

^tems of this liowever, are usually extremely complex: and costly v 
♦ 

An alternate approach to this prc*>lem is to allow the balloon plat- 
form to swing freely from its suspension point and then enplcy some 
method to determine its attitude (orientation) , The orientation par- 
ameters for the platform are determined ty fitting tJie reSult*5 obtained 
from a mathanatical model (vdiich simulates che balloon system) to those 

t 

results cA)tained from the platform's sensors (i.e. , gyroscopes). 

Several nimericcil parameter estimation methods have been develcped 
to determine the attitude (orientation) parameters. Ilie problem is 
normally solved by esiploying an optimization process \duch minimizes 
the error between predicted and known output results. In the case of 
balloon research platforms the cptimization problem involves the mini- 
mization of the sum of the squares of the differences between the angu- 
Icu: velocity ooiponents obtained from the rate gyrosoopeS and those 


4 


# 


predicted fran the nathanatical model . However, with this approach, 
the problem of determining the qptimal decision variables, (i,e*r •iJii” 
tial condition parameters) can require considerable ocnputer running 
time. 

A premising, new i^roach for determiniiig attitude of ballocn re- 
search platforms involves observer state space reconstruction. For 
totally observable systems, state estimators can be constructed. The 
state estimator, \diich is driven by all pleuit inputs and outputs, can be 
used to determine the system state . 

Observer systems are state estimators constructed such that the 
error in the estimated state decays to zero over a finite time interval. 
Ey subtractii^ the plant model frem the chserver model, the ^rror model 
for the reconstructed state can be determined. This model consists of 
a system of homogeneous, first order differential equations. The eigen- 
values of the resulting eigenvalue problem can be chosen such that the 
error decays to zero in a small interval of time. In this case, then, 
the estimates response approaches the actual states exponentially. 

1.3 Object of Thesis 

The two main objectives of this thesis are given as follows; 

1. Develop an observer model for predicting the orientation of 
balloon home research platforms. 

2. Qiplcy this observer model in conjunction with actual data ob- 
tained from NASA'S LfiCATE mission in order to determine the 
platform orientation as a function of time. 

In order to achieve the above objectives it will be necessary to 
first develop a general three dimensional mathematical mcxiel for simu- 
lating the motion of the balloon platform. This will be discussed n^. 
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CHAPTER n 

DEWEDdWEJW OF BAUtJOON SYSaTM MKIHEMATICAL MODEL 
2.1 Ide^dizlltions of System 

The general balloon system to be studied in this report Is shot'n in 
Pig. 1,1-1. The actual motion Of this i^stem is very ocnplex and in- 
volves various iypes of oscillations including bounce (vertical oscil- 
lation) , pendulations (inplane notion) and spin (rotation) . In general, 
it is neoess 2 uy to first idealize this systan before developing the 
mathematical model. For purposes of this study, the following idealiza- 
tions will be made: 

1. The mass of the balloon, subsystems and lnt<aroorj« 2 Cting sub- 
systems, will be ''l\snp©l” at the locations shown in Fig. 2.1-1, 

2. The balloon will be treated as an "equivalent” rigid body. 

3. The altitude of the balloon static equiXibrivin position (float 
atlitude) will be assumed to be a constant during the entire 
period of observation; i.e, , changes in this altitude due to 
losses or changes in the properties of helium will be neglected. 

4. Ihe interconnecting cables will be considered to be; inflexible. 
The above idealizations were ^lied to the general, balloon system 

shewn in Fig. 1.1-1. The resulting idealized systom is shown in Fig. 

2 . 1 - 1 . 

There are two alternate approaches which can be followed for pur- 
poses of modeling the balloon system; these are svznnarized below. 

1. The mathematical model for the entire ballocxv system (Pig. 

2.1-1) can be developed. The major disadvantage of this ap- 
proach is that it requires knowledge of the aerodynamic forces 
acting ui the ballocn itself. Moreover, with this model, a 





nunber of generalized ax>rdinates are needed to specify 
the configuration of the system. 

2. An alternate approach is to develop the mathematical model for 
predicting the motion of subsystems one and two (Pig. 2,1-2) . 

This model does not Include the aerodynamic forces acting on 
the ballcon; hcx^ever, it does require information on the motion 
of the i”:adar reflector support point 0 (Pig. 2,1-2). The 
advantage of this model is that the nuniDer of degrees of freedom 
is decreased, and the aerofynamic force effects are autcmatically 
included if the motion of the support point is knewn. This is 
the model vhich will be used for carryiiig out the researdi dis- 
cussed in this thesis. 

2,2 Generalized Coordinates 

The generalized coordinates for a given system are those coordin- 
ates vhich are enplqyed to specify the configuration of the system at 
any instant of time. In any mechanical system there will be as many 
generalized cxxjrdinates eis there are degrees of freedom. In the case of 
the idealized lumped system shown in Fig. 2.1-2, six generalized coordin- 
ates are required to specify the balloon configuration. These are com- 
prised of six Euler angles vhich specify the orientation of the two 
subsystems. The three translational coordinates located at the radar 
reflector s^jpport point 0 are not considered to be generalized coordin- 
ates, since these are )ciown (presexibed) fixm data obtained fron the 
radar tracking installation. 

In general, the Euler angles give the orientation of the body co- 

1 1 f 

ordinate axes (X . ) relative to a fixed coordinate system (X, ) . A 
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series of three rotaticMis about the bod/ axis is sufficient to allcw the 
bod/ axes to attain any desired orientation. 

Several sets of Euler angles are possible for fixing the orientation 
of the subsystems. One set of Euler angles (Set III) was enployed in 
this work and is shown in Pig. 2.2-1, The sequence of the three rota- 
tions vdiich define this set is described below. 

a. a posiave rotaUon ^ about the X 3 axis xesultleg in 'the 
body system, 

b. a positive rotation 6 .ibout the axis resulting in the X^' 
body system, and 

o. a positive rotation i), about the Xj' axis resulting in the ’ 
body systesn. 

Ihe transformtion equation for the above sequence of rotatiois is 
given as follows; i.e. . 


— .• • • — 

X - AX 

vdiere 


(c(ii;)c((zJ)-s(0)s(0)) (c(t|))s(<zJ)+s(e)c({z^)s(4»)) (-s(qi)c(0)) 

{■<l{Q)Si9i)) (G(e)c(0)) (S(0)) 

(S(ilj)C(jzi)+S(e)S({zJ)c (ilJ)) (S(jzOS(4j)-S( 6)C (J^)c(i|j)) (c(0)c(qj) ) 

X denotes the fixed system axes, and 
X' ' ' dianptes the fixed body system axes. 



2.3 Lagrange’s Equation 

The mathematical model for simulating the motion of ■the research 
platform will be developed by enplqying Lagrange’s equaticai. The general 
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form of this equation is giv«i as follows; i.e. , 


d 3L 
dt 


3q^ 


1 ^= 0 . 


(i “ If » • # f n) 


where 


(2.3-1) 


L - T-V « Lagrangian, 

T « kinetic energy of the system, 

V * potential energy of the system, 

= generalized coordinates, 
q. = generalized velocities, 

mL 

n = the nuitiaer of generalized coordinates, and 
= the nonconservative generalized forces. 

For purposes of this work the friction at the support points 0 and 
1 will be neglected. In addition, the aerodynamic drag forces acting on 
sxibsystans 1 and 2 will also be neglected. Hence, the generalized forces 
are equal to zero and Eq. 2.3-1 reduces to 


d 3L 3L 
dt ^ ” 3g . 

9% 


= 0 . 


(2.3-2) 


2.4 Kinematics 

In order to c±>tain the kinetic energy of the system, it is first 
necessary to develop ‘the kinematic expressions for the velocity (angular 
and linear) of the subsystems. By enploying the Euler angles, the 
angular velocities for subsystem (1) are given as follows; i.e. , 

= (-i^ c(ej^) + 0;^ cw^)- 


„ (1) . 


' 2 ' ' = S(9j^) + and 


(2.4-1) 


n ) • 

^3 = ({zJiC(e^)c(4^) + s(ijjj^). 


. i3 


where 


ORlQiNAt PAQE IS 
OF '^R QUALITY 

« ocirponent of the emgular velocity of sub^stem (1) 
along the i th body axis (i » 1^2,3) , and 
®1'^1'^1 * angles of rotation for subsystem (1) . 


Similarly, the angular velocities for subitem (2) are; 


Wi^^^ = (izl2 S#2)C(e2) +02 0(1^2)), 

« (^2 S(«2^ + ^ 2 * (2.4-2) 

* (^2 C ( V ‘^'( 4 ^ 2 ^ + «2 ^^^^ 2 ^' 

vdiere 

( 2 ) 

- oarponent of the angular velocity of subsystem (2) 
eilong the i th body axis (i - i#2,3) , and 
= Euler angles of rotation for subsystem (2) . 


By enploying small angle approximations (i.e. , S(6) = .0, C(0) = 1) , 
and by neglecting second order terms in 0 and , Etf • (2 . 4-1) and (2.4-2) 
can be written as follows: 


and 


’1 

W 

W 

'^3 




W ^2) 

W <2) 

W <2) 
3 



(2.4-3) 


(2.4-4) 
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Ihe translatlowd motion of the support point 0 (Fig, 2,1-2) is 
referred to m axis which is ^ixed in space. The translational vel- 
ocity expression for this point is given as follows; i.e., 

(2.4-5) 


V3<0>.x; 


vdiere 


V. 


( 0 ) ^ 


absolute velocity oonponents of support point 0 

(i « If 2, 3) along the X. ' " body axis of subsystem (2) 

3* 


The velocity expressions for support point (1) (Fig. 2.1-2) are 

». 

given as follows; i.e., 


V, 


V. 


( 1 ) = 
( 1 ) = 


(Xi - , 

(X2 + 

(ij), 


(X 2 + 


(2.4-6) 


v^ere 


= absolute velocity oonponents of point (1) along the 
boc^ axis of subsystem (2 ) , and 
= distance between point (0) and point (1) . 

The velocity expressions for point (2) are given as follows; i.e. , 


V, 


V, 


V, 


(2) _ 


( 2 ) _ 
( 2 ) .. 


~ (^x "" ^x'^ ^ 

(X 2 + + r292 ) , and 


(2.4-7) 


v^ere 


(2) 

V. ' ' absoltite velocity conponents of point (2) along the 


bod/ axis of subsystem (2) , and 
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■ distance between point (i-1) and point (i) . 


The above equations were developed by errploying small angle aqpprojc- 
imations (i.e., S(0) « 6, C(0) * 1), and neglecting second order terms 
in 0 and Moreover, it also assumes that the nature of the ring and 
clevis s\jpport at point (1) is such that the difference between the spin 
angles 0^ and ^2 small. 

2.5 System Lacrcangian 

The general kinetic expression for subsystems (1) and (2) is given 
as follows: 


where 






(i=l,2) 


(2.5-1) 


T s= kinetic energy of subsystem (i) , 
m. * mass of subsystem (i) , 

■lo 

\l' ^i2' ^i3 " of inertia of subsystems (ij^ along 

the (2) body axis, and 

vfj s= oonponents of the angular velocity of subsystem (i) 
cdong the (2) body axis. 


The totcil kinetic energy T of the balloon system is obtained by 

1 

sunning Eg. (2.5-1) and substituting from Bqs, (2.4-3), (2.4-4), (2.4-6) 
and (2.4-7); this gives: 


T = T^^^ + T^^^ 

+ » 3 >^> 

+<VW’^2«2>^ + (Xj)^) + ^ 


(2.5-2) 
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For purposes of developing Eg. (2.5-2) the moments of inertia 
^22 neglected ’’and and m 2 were treated as point masses. 

The systefii potential energy is due to the presence of the conser- 
vative gravitational forces and is given as follows: 

V » (2.5-3) 

vAiere 

* -m 2 g(r 3 ^C(ej^)G(ijJj^) + r 2 C(e 2 ) 0(1^2) ) • 

Ihe system Lagrangian (L ) , which is defined in Fq. (2.3-1 ) , is 
cdotained by subtracting the total potential energy (Eg. (2.5-3)) from 
the total kinetic energy (Eg. (2.5-2)) and is given as follows: 



+ mi griC(ei)c(i()i) 

+ m^ g (riC ( 61 ) c (il^i) + r 2 C (62) C (4/3) ) • (2 . 5-4) 

2.6 ^stem Math Model 

The equaticais for the motion of the ballocn platform are obtained 
by substituting the Lag’rangian fran Eg. (2.5-4) into Eg. (2.3-2). The 
resulting equations are given below. 
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( 2 . 6 - 1 ) 


( 2 . 6 - 2 ) 


\^iere 


« Or and 

(2.6-3) 

•• 

^2 « Or 


= (m^ + 



"^1 “ ”^V2' 
m22 « ^2^ 


’'u “ 

kjj = iti2 grj. 



^2 ~ V2' 

= acceleration ccm«nent o£ point (0) along the e/^) bo*r 
axis, and 

( 2 ) * 

a 2 = acceleration cxatfxDnent of point (0) etLong the e^^' ' bo(fy 
axis. 


Bgs. (2.6-1) and (2.6-2) were developed assuming small displaoements; 
l.e. , C(6j^) “ C(4^j^) ~ Ir “ ®j_r S(4^j^) ~ 

Eg. (2.6-3) yields that 96 ^ is a constant. In this study the pre- 
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cision of this :npdel was inproved by eni>loy*4ig the ixansfornation ogua- 
a™ for the angular velocity . fta atated earUer, it will be aa- 
Burned that » p< 2 * ^^us. Eg# (2a6-3) will be replaced with the follcw^ 
ings 

« ^2 * ^3^^^ (2.6-4) 

v^iere 

W 3 < 2 ) is the angular velocity obtained from the rate gyro 

( 2 ) 

mounted along the e^ body axis. 

Bgs. (2.6-1) and (2.6-2) can be written in matrix form as follows? 


where 


Ei-Chi 


hi 


^2 = 


M 


K ~ 


i = Vi 

(i=l,2) 


(2.6-5) 



(i“1.2) 


(2.6-6) 






«2 

t 









^2 

t 




^1 

"^1 

”'l2 

"^2 

t 



P^ll 

0 ” 





^22 



1 ? 
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W,(2) 


F 



« 

B 0 


SB 

2 


the ooitfionent of the angular Yelocity of subeystem 
( 2 ) 

(2) along the bo^ axis, 


W, 


(2) 





the cotpon^t of the angular velocity of subsystem 
(2) adong the © 2 ^^^ body axis, 


C » (0 1), 



, and 


nie nmerical values fa the n.y and ky coefficients were oorputed 

4 

fcy employing the properties of the balloon system vdiich are given in 
Table 4.1-1. These resulting values are given in Table 4.1-2, 

2.7 State Variable Form of Math Model 

Egs. (2.6-5) and (2.6-6) can be expressed din general state variable 
form as follcws: 

L q = Nq + R u, and (2.7-1) 


y - C q 


(2.7-2) 


vyhere 





u = cxjntrol variable, 



= the cxinpcsnent of platform angiilar velocity along the 
astern (2) body axis, and 

C = (0 0 0 1) . 

Bqs. (2.7-1) and (2.7-2) csan be further expressed as follows; 


q = Aq + Bu , and 


(2.7-3) 
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v^here 


y « C q 
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(2.7-4) 


A * l“^ N 


0 

0 


a. 


0 

0 


31 '‘32 

‘41 .^42 


1 

0 

0 


0 0 


B s l"^ R = 


^ elements of the A and B mtrioes are given below? i.e. , 


®31 “ 


_ -™22 ^11 


^32 " 


_ "^12 ^22 


_ "^1 ^3.1 
41 j ' 

“"h.1 ^22 


^42 = 


^3 = 


"^2^1 " "'l2^2 




v4iere 


j 1^22 - ii^2i"\2' 


Ihe numerical values for the a^j and coefficients were c c np uted 
by employing the values of and m^^ vAiich are given in Table 4.1-2. 
Ihiese resulting values are given in Tables 4.1-3 and 4.1-4. 
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DEVELOPMENT CP OBSER/ER SYSTE*! MRTHEMATICAL MODELS 
3.1 Concept of Ctoservability 

Ihe observability of a system iirplies the deterrninability of the 
system state frcm an observation of the output over a firiite time in- 
terval starting from the liistant at vdiich the state is desired' 1 It 
is assumed that the system inputs, outputs and mathematical model cire 
kncwn. 

For purposes of referencing the work in this chapter, the state 
variable form of the balloon's model (Bqs. (2.7-3) and 2.7-4)) will 
be rewritten below; i.e. , 

q = A q + B u, (3.1-1) 

q (V 

y = Cq. (3.1-2) 

where 

q is the n th order state vector, 
is the unknown initial state vector, 
u is the single (scalar) input, 
y is the single (scalar) output, 
n is the order of the system, 

A is a n X n matrix, 

B is a n X 1 matrix, and 
C is a 1 X n matrix. 
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Bq. (3.1-1) determines the plant (Ramies and indicates how the 
input (or control) u affects the state vector q. The matrix A char- 
acterizes the plant dynamics when u is not present; this is the so- 
called free-response case. The matrix B determines how the plant re- 
sponse is affected by the input (or control) vector u. Eg. (3.1-2) 
defines the relationship between state vector q and the syston output y. 

For system ebservability, the basic question is as follows: "Is 

it possible to identify the initial state q^ by observing the output 
(y) over a finite time interval?" Precise definitions of system ob- 
servability are given eis follows: 

1. Definition 3.1a. A state q^, i.e. q(tj^) of a system is said 
to be observable at time t^, if knowledge of the input u(t) 
and output y(t) over a finite time t^ < t < t^,. ooTpletely de- 
termines the state q^. Otherwise, the state is said to be 

unobservable at t^, 
o 

2. Definition 3.1b. If all systen states q(t) are observable, 
then the system is said to be ocrpletely observable or just 
observable. 

3. Definiticm 3.1c. If the state q^ is observable and if the 
knowledge of the input and the output over an arbitrarily small 
interval of time suffices to determine (independent of t^) , 
then the state is said to be totally daservable. 

4. Definition 3. Id. If all the states q(t) are totally c±>servable, 

if, 5) 

then the system is said to be totally observable,' . 

The necessary ocxidition for observability of the balloai systen is 
given in the following section. 
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3.2 Observ^ility of Balloon Systems 

The follofr/ing theorem (Ref. 3) can be enplqyed 'to determine the 
observability for general time - invariant ^sterns. 

Theorem 3.2a. The time-invariant system described by Eqs. (3.1-1) 
and (3.1-2) is totally cdaservable if and only if the ooTfosite matrix M 
has rank n, v^here 


C 

CA 

.C(A^) 

c(a"-1) 


and 


n, C and A are defined in Eqs. (3.1-1) and 3.1-2). 

The observability of the balloon system described by Eqs. (2.7-3) 
and (2.7-4) can new be determined by showing that the oonposite matrix 

2 3 

M has rank 4. For this purpose, the matrices C, CA, C(A ) and C(A ) are 
given as follows: 


C = (0 0 0 1.) , 

CA = (a^j^ a^2 0 0)/ 

C(A^) =00 

C(A^) = (( 34133 ^^ ' ^42^41^ ^^32^41 ^ * 

The specific form of the oonposite matrix M for the beilloon system 


model is given as follows, i.e. 


M = 


CA 

C(A^) 

C(A^) 




; i.e. 


l|r 
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M = 


‘41 

0 




42 

0 


^^ 32 ^ 41 ‘^^ 2 ^ 


0 

0 

^41 

0 


1 

0 

^42 

0 


In order to show that the rank of M = four, it is necessary to 
prove that the determinant is non zero. It can be shown that the de- 
terminant for M is non zero if the following expression is non zero, 
i.e., if 

Thus, the balloon system is octrpletely dsservable since the ocxnposite 
matrix M has rank = 4. 

3.3 Cfcservability of fealloon System With Output Bias 

Previous studies have been conducted for the LACATE :^stem in order 
to determine the natxare of the balloon's platform motion The time 
history of t!,ie platform pendiiLation angles (6{t) and ijj(t) ) was deter- 
mined by integrating the output fran the rate gyros. Ihe stu^ indi- 
cated that the platform motion consisted of small oscillations sii^Jerim- 
posed on a line with (nearly) constant slope. These results suggested 
that the gyroscopes contain a constant bias error. 
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In order to take into account* the error in the output (y) due to 
bias in the gyrosoopes, the matrices in Bqs. (3.1-1) and (3.1-2) are 
defined as follows; i.e., 



u = a, 
y = w‘2> + 

( 2 ) 

w = the ocfiponent of the angular velocity of the platform 
along the specified body axis, 

qg = corresponding gyroscc^ bias, and 

C = (0 0 0 1 1) . 

The cbservability of tlie system model described by Eqs. . (3.1-1) 
and (3.1-2) can be ascertained by showing that the corresponding ocin- 
posite matrix M has rank 5. It can be shown that the ccirposite matrix 
M for this system is given as follows; i.e.. 
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c 

CA 

C(A^) 

C(A^) 

C(A^) 
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M = 


0 

^41 

0 


0 


(^Al +a^^a 

ill 


*41“31 “42“4X' 
0 


0 

0 

^41 

0 


a 


1 

0 

42 

0 


1 

0 

0 

0 


^^41^3l‘^^42^41^ ^®41®32'^®42^ ° 


It can be shown that the detenninant of the above matrix is non-zero 
if the resulting expression is non zero; i.e. # if 

nijm2 ^ 

Uius, the system described ty Bqs. (3.1-1) and (3.1-2) is oonpletely 
observable since the ccmposite matrix M has rank = 5. 

3.4 Development of Full Order Identity Observer for Balloon System 
Without Bias 

An n th order identity observer (or asynpotic-state estimator) can 
be constructed for the ocnpleteiy observable n th order plant described 
by Bqs. (3.1-1) and (3.1-2). The observer is described by the following 
equations; i.e. / 

# 

Z = FZ + B u. + G y, and 


(3.4-1) 


v*iere 


Z is an n th order estdmate of the state vector q, 

Zq is an estimate of the unkrKJwn initial state vector <3^, 

G is an n X 1 matrix, 

F = (A - GC) , and 

y, u, B and C are as described in Egs. (3.1-1) and (3.1-2). 

An inspection of Bq. (3.4-1) reveals that the state estimators 
response (Z) will be determined from a consideration of the observer 
dynardcs, eictemal inputs and plant outputs, The observer dynamics are 
controlled by the P matrix. Ihe external input u contributes to the 
state estimator's response via the oontrol matrix B. From Bq. (3.1-1) , 
it can be seen that this control matrix B and the external inputs u are 
identical for both the plant and observer. 

For accurate state space reconstruction the plait output y must be 
fed into the observer model. By coupling the plant output y to the ob- 
server via the G matrix, the observer becomes a closed loop estimator. 
This is illustrated in Pig. (3.4-1) . 

The P matrix in Bq. (3.4-1) is constructed such that the difference 
bebyeen the output of the observer mo^l and the plant model is zero 
over some finite time interval. This error (E) is given cis fellows; i.e. 

E = Z - q. (3.4-2) 

Subtracting Bq. (3.1-1) from (3.4-1) yields the following: 

E - FE. (3.4-3) 


The solution of Bq. (3.4-3) yields: 




Fig. (3.^-1) Block Diagram For Asymptotic State iSstlma tor 
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E s E 
o 


(3.4-4) 




®o ' - %■ 


It is clear from Eq. (3.4-4) that the error, E, decays exponeitially 
to zero if G is chosen such that all of the eigenvalues of P are nega- 
tive or have negative real parts. Also, these eigenvalues must be more 
negative than the eigenvalues of A to insure accurate response. 

ihe detailed form of the F matrix in Bq. (3.4-3) car- be obtained by 
substituting the form of the A and C matrices vhich are defined in Eq. 
(2.7-3) and (2.7-4). Uiis yields the following; i.e. , 


P — (A — GQ) , i.e. , 

" 0 0 1 -g^“ 

0 0 0 ( 1 - 92 ) 

®31 ^32 ° "^3 

.41 .42 0 -^4 

where 



Hie form of the solution for E in Bq. (3.4-3) is given as follows 
E = X e^^. 
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Substitution of this into Eg, (3*4-3) yields the following eigenvalue 
problem: 


(P - XI)X « 0 

where 

I « the identity matrix, and 
P « matrix defined in Eg. (3.4-1) , 


Tlie necessciry and sufficient cxDnditicMi for determining the eigen- 


values of the matrix P is 

that 

> 

1 F-XI 

: 1? 

i.e. , 



-X 

0 

1 

"^1 


0 

-X 

0 

(i-gj^ 


^31 

^32 

-X 

"^3 


^41 

^42 

0 

r< 

1 

Ihe characteristic polynomial 


given belcw? i.e., 

+ <(- 94 ® 31 '^ 3 ® 41 >^>+<'® 32 ' 
+ g2®32®4l'®2®42®31* “ ® 


* 0. (3.4-5) 

obtained by expanding Eg. (3.4-5) is 

'4l'^^42^31 

(3.4-6) 


Critical damping of the error E is obtained by determining the G 
matrix such that all of the eigenvalues are negative and equal. Ihis 
yieldi? the following values for the G matrix: i.e. , 
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■ «-6^Si-»42^2“42>/(-“41»' 

92 * ((■^^■®32®41''’®42“31*/*““32“41’*'®42‘'3i'*' 

93 * , and 

94 ® (""4A.) » 


3.5 Develogment of Identity Ctoserver Systems for Balloon System With 
Biets 

In the case of the flf ih order identity observer system the matrices 
in Eg. (3.4-1) have the following form; i.e, , 


0 

0 

1 

•^1 

“^1 

0 

0 

0 

( x - g 2 ^ 

"^2 

^31 

^32 

0 

“^3 

"^3 

^41 

^42 

0 

-^4 

•^4 

0 

0 

0 

■^5 

"^5 


G = 


and y, u, B and c are as described in Eqs. (3.1-1) and (3.1-2). 

The eigenvalues for the P matrix corresponding to the biaus model 
are obtained fron the foilwing oonditicai; i.e., 


^1 

^2 

% 

^4 

% 
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-X 

0 

1 

-^1 

"^1 

0 

-X 

0 

( 1 - 92 > 

■^2 

®31 

®32 

-X 

^3 

■^3 

®41 

*41 

0 

(-^-54) 

■^4 

0 

0 

0 

-95 

(- X - gg ) 


Olie characteristic polynanial of Bq. (3.5-i) is given as follows? 

i*e* f 

(X®)+'(g4+gg)X^+( <gia4i+?2»42''‘3l'®42)^^> 
•"<<53»4I-^4®31-^5®31-^5®42>^^> 

+ ( (-g2*41®3l'*«2®41^32+“42®3r*41®32> 

’*'^%^42^3l“^5®41^32^ ” (3.5-2) 

The final form of the G matrix elements (obtained from the ocxidi- 
tion for critical damping) are given as follows? 

gi “ ^^^°^^“^2^42'^^31‘^^42^'^^41^ ' 

gg * USX *a42®3i‘^®41®32^/^”^42'*3l’*’^41®32^^ ' 

93 = (KOX^^g^aj^-tgj ( 833 + 3 ^ 2 ) )/a^ 3 >, 

* 

?4 = (“*51. -gg) , and 

95 = ( (-X )/a4233l“'a3^-a^j^a22) ) • 
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CHAPTER IV 

RESUnrS AND OONaUSIONS 

» 

4.1 Data for LACATE Mission 

Figure (1,1-1) illustrates the actual LACATE balloon syjstem and 
Figure (2.1-2) illustrates the corresponding idealized system used in 
this stu(ty. The values for the various lengths and inasses of the ide- 
alized sysicem are given in Table (4.1-1) . 

Hie numerical values for the elements of the M and K matrices of 
Bq. (2.6-5) were coiputed based cai the data given above. The resulting 
values are presented in Table (4.1-2) . The numerical values for the A 
and B matrices of Eg. (3.1-1) are given in Tables (4.1-3) and (4.1-4). 

The eigenvalue problem for the balloon system was solved analyti- 

2 

cally. The solution for the eigenvalues (Q ) j aixi corresponding eigen- 
vectors are presented in Table (4.1-5) . The values of Qj represent the 
natural frequencies of the system. The modal sh^Je functions and per- 
iods corresponding to each natural frequency are shewn in Table (4.1-6) . 

Results for the balloon observer system were obtained by employing 
two separate time inter'/als. The equations for predicting the boefy axis 
acceleration coiponents from sensor outputs are presented in Appendix A. 
Plots of the body axis acceleration components over the two time inter- 
vals are given in Figures (4.1-1) through (4.1-6). 

4.2 Results From Simulation Stuc^ 

Bq. (3.1-1) was enployed to simulate the balloon platform angular 
velocity and angular displacement. The inputs and outputs employed with 
the simulated system were of the same order of magnitude as predicted 





II 


TABLE 4.1-1 

Idealized LACATE System Properties 

(distance from point 0 to mass = 75 ft. 

T 2 (distance from mass m^^ to m 2 ) * 15 f;t. 
m^ (limped mass) = 135 

”"2 = 375 ib^ 


TABLE 4.1-2 

Coefficients of m and k Matrices 

m^^ = 89156 (Ibf • • ft) 

m^2 = 13111 (Ibf • s^ • ft) 

1022^ = 13111 (Ibf • • ft) 

IH22 = 2622 (Ibf • s^ • ft) 

= 38278 (Ibf • ft) 
k22 « 5629.4 (Ibf • ft) 
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TABIE 4>1"3 

Coefficdents of A Matrix 


* -1.621 (s"^) 
= 1.1923 (s“^) 
=8.107 (s”^) 

= 8.107 (s"^) 


TABLE 4.1-4 


Coefficiaits of B Matrix 


= -0.01334 (ft”^) 


= 0.0 


= 0.01334 (ft“^) 
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Pig. (^.1-3) Plot of Balloon's Translational Acceleration (a^) (0=tt«500) 




Fig. (^.1-4) Plot of Balloon's 









by the ac±ual lACATE flight data. For cxnparison pv»:poses, all of the 
3irnu3.ation runs eni>loyed identical initial conditions and eigenvalues; 
the magnitude of the latter was equal to ^0.6. 

The angular velocity (0) and angular displacement (0) predicrtied by 
the fourth order observer system are shown in Pigs. (4,2-1) - (4.2-4) 
for the case vtei bias is not present in the output. The free response 
case is iXl';istrated in Figs. (4.2-1) and (4.2-2) , v4iile in Figures 
(4.2-3) and (4.2-4) , results are given for the case vdien a ranp input 
was cnplcyed, 

A ■ 

Figures (4.2-5) - (4.2-13) illustrate the free response of the 
fourth order observer systan for the case vdien bias is present in the 
output. The plant output for the case shown in Figs. (4.2-5) and (4.2-6) 
cx>ntalns a constant bias, v4ule a linear bias was used to cfctain the re- 
sults shown in Figs. (4.2-7) and (4.2-8) , 

Figures (4.2-9) - (4.2-13) give the free response results for the 
fourth order observer system when a high frequency bias = 0.005 
cos(l3t)) is present in the plant output. The results for the angular 
displacement from the fourth order observer system are presented in 
Figs. (4.2-9) and (4.2-10). Figs. (4,2-11) and* (4.2-12) present the 
results for the angular velocity predicted by the fourth order observer 
model. Fig. (4.2-13) presents the results for the plant output (i.e., 
y = q4 + 3). 

Figures (4.2-14) - (4,2-21) display the free response results from 
the fifth order observer system for the case v^ien bias in the output is 
present. In the case of Figs. (4.2-14) and (4.2-15), the output con- 
tained a constant bias, vrtiile a linear bias was used to obtcdn the re- 
sults shown in Figs. (4.2-16) and (4.2-17). 


^5 


Figures (4.2-18) - (4.2-21) show the results frcm the fifth or^ 
observer system When the plmt outp(it contiiins a high frequency bias 
(qg « 0.005 cos (13t)). Pigs. (4+2-18) and (4.2-19) present results 
for the angular displacement predicted by the fifth order <A>server j^s- 
tem. The angular velocity' predicted by this system is i^icMn in Figs. 
(4.2-20) and (4.2-21). The results for the plant output (i.e., y * 

+ qg) are presented in Pig, (4.2-13). 
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ilts For Aiig\ilar Displacement 

er, w=0. >i=-0.6. Bias=0,00^cosfl^t')’i 



Fig. (4.2-19) Simulation Eesiilts For Angulax Displacement 

(5th Order Observer, u=0, /\=-0.6, Bias=0,005cos(l3t)) 
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Fig. (4.2-20) Simulation Results For Angular Velocity 

(5th Order Observer, u=0, A=-0i6, Bias=0,005cos(l3t)) 
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4.3 Discussion of Simulation Study Results 

Hie results in Figs, (4. 2-1) -(4. 2-4) verify that, after a finite 
time interval, the fourth order observer system will reconstruct the 
plant state exactly if all plant inputs and outputs are known. The 
use of repeated eigenvalues (i.e, , - - 0.6) resulted in an error free 

response after a period of 20 seconds. 

The fourth order observer results in Figs. (4. 2-5) -(4. 2-13) were 
obtained for tlie case vhen the plant output contained a bias (3(t)) . 

From Figs. (4.2-5)-{4.2-13) , it is seen that the reconstructed states 
contained an error resulting fron the bias presoit in the plant output. 

TVjo methods can be used to determine the steacty state error in this 
reconstructed state due to output bias: 

1. Revise Eg. (3.4-3) to include the error due to bias and then 
solve this equation. 

2. Ccnpare the results for the reoonstructed state to the actual 
state directly by employing the curves in Figs. (4. 2-5) -(4. 2-13) . 

For purposes of this work the latter metlKxi will be enployed. 

For conparing the reconstructed stat ,'s to the actual state, a relative 


error, ER^(t^) will be used. This is defined as follows; i.e.. 


ERj^(tj) = 


EA^(tj) 

q]7t^ 


(4.3-1) 


where 


EAi(tj) = magnitude of the error in the state at t = tj, and 


qi(t.) = maximum value of the actual state variable at t = t.. 
The curves in Pigs. (4.2-5) and (4.2-6) indicate that the error in 
the reconstructed state is a constant vhen the plant output contains 
a constant bias (i.e., 3 = 0.001) . In particular, for this case, a 5% 


4 
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relative error (ER ) in the output resulted in errors (EEtg and ER^) of 

** • 

17% and 100% in and 9 ^ respectively. 

—6 

When a linear bias (i.e., 3 = 4.0x10 t) is present in the plant 
output, a linear error resulted in the reconstructed state (see Figs. 
(4.2-7) ani (4.2-8)). At t = 43s, ER^ = 6%, ER 2 = 23% ani ER^ = 133%. 

Ihe results shown in Pigs. (4. 2-9) -(4. 2-13) indicate that vrtien the 
plant output contained a high frequency bias (i.e., $ - 0.005 cos 
(13t) ) , the angular displacement (9) and the angular velocity (9) pre- 
dicted by the fourth order observer oontair^ a sinusoidal error. The 
error in the reconstructed state at peak arrplitudes cure ER^ = 167%, ‘ 

ER 2 = 9« and ER 4 =. 33%. 

The results obtained from the fifth c^der c^server are shewn in 
Figs. (4. 2-14)- (4.4-21) . Figs. (4. 2-14) -(4. 2-16) indicate that, for the 

4 

case when the plant output contains a constant bias, the actual state will 

be reconstructed exactly after a finite interval of time. 

Figs. (4.2-16) and (4.2-17) show that viien the plant output contained 

a linear bias (i.e. , $ = - 1 x 10 ^t) , a constant error resulted in the 

reconstructed state. At t = 43s., the resulting errors are given as 

follows; i.e., ER^ = 153%, ER 2 = 64%_ and ER^ = 100%. 

Pigs. (4. 2-18) “(4. 2-21) show that when the plant output contains a 

high’ frequency bias (i.e., 3 = 0.005 cos (13t)) the state predicted by 

the fifth order observer contains a sinusoidal error. The errors in the 

reconstructed state at peak aitplitudes are ER^ = 167%, ER 2 = 6% and 

ER. = 35%. 

4 

In general, the results of the simulation study shew that the form 
of the error in the reconstructed state will depend both on the order 
of the observer model and the type of bias (i.e., constant, linear, etc.) 
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present in the plant output. VJhen an nth order bias is present in the 
plant output, the fifth order observer will yield a more accurate response 
than the fourth order diserver. For this case, the error in the fifth 
order ctoserver will be of order n-1, while the error in the fourth order 
observer will be of order n. 

For a specific form of bias, the magnitude of the errors in the re- 
constructed state will depend on the elements of the F and G matrices* 

For the particular case when - - 0.6, the magnitude of error in the 
reconstructed state Z 2 was much smaller than that in Z^. Moreover, for. 
the case vhen the bias was of a sinusoidal nature (B = 0.005 cos (13t)) , 
the error in the reconstructed state 7>2 fifth order 

observer model; hcwever, a more accurate response was dbtained in state 

r 

Z. for the fourth order observer. 

4 

i 

4.4 Observer Results for Determining Orientation of Balloon Platform 
The fourth and fifth order observer system models were enployod 
to determine the orientation (€ 2 ) of the balloon platform in the X 2 X 2 
plane for the time interval t = Qs (initial LACATE data recording time) 
to t * 500s. The observer transient error was assumed to decay to zero 
after an elapsed time period of 250s. At this time the initial condition 
for integrating the gyrosoc^ was set eqiaal to the angular displacement 

(Z 2 ) predicted by the observer. Fig. (4.4-1) illustrates the output (y) 

1 1 1 

datained from the gyroscope with sensing axis along the e^ bo(^ axis. 

Pigs. (4,4-2)-(4.4-5) give the angular velocity (Z^) and angular 
displacement iz^) predicted hy the fourth order observer nodel for the 
case vdien = -0.5. The free response case is illustrated in Figs. (4.4-2) 
and (4.4-3), while Pigs. (4.4-4) and (4.4-5) give the response when the 
wind acceleration is included. 
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Pigs. {4.4-6)-(4.4-13) illustrate tJie free response of the fifth 
order observer i^stem with =-0.2 and - 0.5 respectively. Figs. (4.4-6) 
and (4.4-10) present the results for the angular velocdty while Figs. 
(4.4-7) f (4.4-8), (4.4-11), and (4.4-12) present the results for the 
angular di^lacement. Pigs. (4.4-9) and (4,4-13) contain the results 
for predicted bias (Zg) , 

Figs. (4. 4-14) “(4. 4-25) present the fifth order observers response 
with wind acceleration for the case v4ien = -0.2, -0.5, and -0.7 
respectively. Pigs. (4.4-14), (4.4-18), and (4.4-22) present the results 
for the angular velocity vrtiile Figs. (4.4-15), (4.4-16), (4.4-19), 
(4.4-20), (4.4-23) and (4.4-24) present the angular displacement. Ihe 
bias predicted from the fifth order observer for the respective cases 
is presented in Figs. (4.4-17), (4.4-21), and (4.4-25). 



Fig. (4.4-1) Actual Output For Gyroscope Counted Along X. Body Axis 
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Fig. (4.4-11) Actual Results For Angular Displacement 
( 5 th Order Observer, u=0, A =- 0 . 5 ) 
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Fig. (4.4-12) Actual Hesults For Angular Displacement 
(5th Order Observer, u=0, A==-0,5) 
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4.5 Discussion of Balloon Observeir Reauita 

Figs, (4,4-2) land (4,4-4) Indicate that the angulai?^leloclty 
results predicted by the fourth order observer differ significantly 
from those obtained from the gyroscope. The actual angular velocity 
of the balloon platform is different than that obtained from either the 
observer or gyroscope. Deficiencies In the ability of the gyroscope 
to reproduce the actual angular velocity are caused by mechanical 

Insensitivity to sudden changes In angular velocity and bias, Errors 

0 " 

in the angular velocity values predicted by the fourth order observer 
result from the form of the elements of the F and G matrices (see 
Sec, 4,3) and the magnitude of bias present In the output of the gyroscope. 
Fig, (4,4-3) shows that the angular displacement predicted by the 
fourth order observer model without wind Input differs significantly 
from that obtained by Integrating the output from the gyroscope (yl). 

The general trend of the angular displacement predicted by this model 
Is simular to that obtained by integrating the output of the gyroscopes. 
However, the output of this model contains errors} e,g,, the maximum 
difference between the two curves Is 0,172°, This Is significant in 

X 

view of the fact that the maximum displacement predicted by the observer 
Is 0,058° while that obtained by integrating the output from the gyro- 
scope Is 0.23°, These errors can be attributed to the fact thati 

1. The effect of the wind acceleration Is neglected, 

2, Bias Is present In the plant output. 

Fig, ( 4 , 4 - 5 ) shows that the angular displacements obtained from 
the fourth order observer Including wind acceleration compare less 
favorably (with the integrated gyroscope output) than those obtained 
from the previous model. The magnitude of the maximum angular displacement 
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predicted by this model is approximately 0.5°» while the maximum 
deviation is of the order 0,6®, 

Figs, (4,4-6), (4,4-10), (4.4-14), (4,4-lB), and (4.4-22) show 
that the angular velocity results obtained from the fifth order observer 
model differ significantly from the results (biased) obtained from the 
gyroscope, This Is true regart\less of whether the effect of wind 
acceleration is Included, However, as the magnitude of the repeated 
eigenvalue Increases, the differences between the two curves decrease. 

Figs, (4,4-7), (4,4-8), (4,4-11), and (4,4-12) show that the angular 
displacements predicted by the fifth order observer model which was 
developed byexcludlng the affect of wind accleratlon, are In good 
agreement with those obtained by integrating the modified output from 
the gyroscope (ylm). The latter’ was obtained by subtracting the bias 
predicted by this model from the actual gyroscope output, These 
figures indicate that the difference between the two curves Increases 
with Increasing time. Moreover, this difference decreases with In- 
creasing magnitude of The maximum displacement wlthXj_“-0.2 and 
-0.5 Is of the order 0.05°{ while the maximum error is 0.08° and 
0,02° respectively. 

Figs. (4.4-15), (4.4-16), (4.4-19), (4.4-20), (4.4-23), and (4.4-24) 
Indicate that the angular displacements predicted by the fifth order 
observer which includes the affect of wind acceleration are in good 
agreement with thv. values obtained by integrating the modified output 
from the gyroscope. The difference between these two curves remains 
constant over the entire observed time period. Moreover, this difference 
decreases with increasing magnitudes of The magnitude of the 
maximum error and displacement for all three cases is approximately 
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0.25° and 0,5° rospRotlvely, 

Figs, (^.^f-5), (^*'.4-15), (4,4-16), (4|4-19)i (4,4-20), (4.4-25), 

<! . 

and (4.4-24) show that the angular displacements predicted "by the fourth 

and fifth order observer models jiompare favorably regardless of the 

' ■ 

value of X. . This Indicates that- the predicted angular displacement 

* ' } 

is unchanged regardless which observer model is used. Thus, the 
fourth order observer model can be used to predict the angular displace- 
ment of the balloon platform. 

Figs. (4.4-3), (4.4-7), (4.4-8), (4.4-11), and (4.4-12) show 

that the angular displacements predicted by the observer models with- 

* • 

out wind input differ significantly from those results predicted by 
the models which Include wind input, -This indicates that knowledge 

ii ■ 

of wind accleratlon Is necessary in order to obtain accurate results 

for the attitude of the platform. 

■' « 

I ■ 

4.6 Qonclusions •' 

This study has shown that, for a completely observable balloon 
system, observer models can be Constructed to accurately determine the 
angular displacement of the observational platform. Any errors in 
the predicted platform state are due mainly to errors in the balloon 
flight data (i ,e, , acceleration and angular velocity data) as opposed 
to deficiency in the observer model. Although the results for the 
angular velocity are In poor agreement with those obtained from the 
gyroscopes, these deviations can be decreased considerably by proper 
choice of the eigenvalues. 

This study has also shoi-m that the angular displacements predlo:ted 
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ijy th® otservir models do not vayy slgnlfloantly with slthoj? the order 
of the model or tho maccnltude of th'd' repeated eigenvalue* However, 

■"'5' “ 

the results do vary slghlfloantly’idepd'adlng on whether the affeot 
of wind aooeleratlon Is Included* 





■flf.PPO** 


100 


APPENDIX A 

Balloon Treunslational AcceXeration Oowtioraaitts 

Xn the case of the lACAH! experiment# the baHoon's position was 
tracked by radar; the tanslatipnal ooqponents were cbtained with respect 
to the earth fixed axis shown in Fig. (A-1) . Ihe oorreeponding body 
axis for the balloon platform system is shown in Fig. (A-2) . Ihe angle 
a measured between these two coordinate systems (Pig. A-2) is given as 


follows; i.e . , 


vhere 


a 


’t 

W 3 dt - tOs(X) 4- (A-1) 

■ o 


W 3 ^ spin ffijqpoimt of angular velocity obtained from gyroscope> 


0 “ magnitude of earth spin, 

(7.2722 X 10“® (rad * s“^)#' 

X *= latitude angle (0.5724 (rad) ) , and 

♦ 

«Q • initial value of a as measured by magnetometer. 

Ihe velocity components of the bedloon were obtained by numerically 
differentiating the (radar tracked) translational components. The velocity 

k 

components of the balloon (V # V ) measured along the balloon’s body axis 

X y 

cire given as follows; i.e., 


where 


- Vj^S(a) + Vy C (ot) , C (a) + S (a) , (A-2) 



= balloon velocity component 
along the 62 bocfy axis , 

Hi 

= balloon velocity ooitponent along the e^ body axis, 

= badloon velocity component along the earth fixed axis, 
s! balloon velocity conponent along the 62 earth fixed axis. 
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an£i Is as daflned in Sq. (A- 1 ) . 

^ btiLloon's translational acceleration conf»nents alongr the body 
axis can be obtained by differentiating Bg. (A-2) with respect to tima. 
The resulting equations are given as following; i«e . , 

32 * V j^C(»)+V S(«)4^(-VjjS(a)+V C(a)), 
vihere 

a^ M translational acceleration acnf)onent along the e2*' ' body axisi 
’ 32 ^ tremslationed acceleration cai|)onent 2dong the e^" ' body aadn, 

» translational acceleration ocnpsnent along the e^^ earth fixed axis^ 

♦ 

« translational acceleration conponent along the 62 earth fixed axisr 


a » W- “ 0S(X ) , and 

J 

V , V and tx are as defined previously. It should be noted thatr in the 

* y 

abo^'o development r the earth's rotaticMial effects are neglected. 
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■1 .000 0 PROGRAM TO OBTAIN BODY ACCEL. ITIOM 

2.000 0 SENSITIVITY ANALYSIS 

‘5.000 OUTPUT ’INPUT N» 

4.000 INPUT N 

5.000 DO 1 I=:1,N . 

6.000 READ(105,2) T,P 

7.000 2 FORMAT (2G) 

8.000 READ(109,2)T,PD 

9.000 READ(104,3) TO, XX, VX, AX 

10.000 5 pormatUg) 

11.000 READ(107,5) T1,YY,VY,AY 

12.000 5 FORMAT (40) 

13.000 A1=AX*C03(P)-*AY#SIN(P)+(-VX*SIN(P)+VYWS(P))*PD 

1 4. 000 A2sAX* (-SIN ( P) ) •♦AY^COS (P) +(-VX*C0S ( P) -VY*SIN(P) ) *PD 

15.000 WBITS(106,4) T,A1,A2 

16.000 4 P0RMAT(3B14.6) 

17.000 1 CONTINUE 

18.000 STOP 

19.000 END 
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2, Balloon Fourth Order Observer 


I .CXX? 0 
2,000 
3,000 

4.000 

5.000 9 
6.000 
7.000 
8.000 
9.000 

10.000 

11.000 

12.000 

15.000 

14.000 

15.000 

16.000 

17.000 

18.000 

19.000 

20.000 

21.000 

22.000 

23.000 

24.000 

25.000 
26.000 1 
27.000 
28.000 
29.000 
30.000 
31.000 
32.000 
33.000 

34.000 

35.000 4 

36.000 5 

37.000 

38.000 

39.000 2 

40.000 
41.000 
42.000 
43.000 

44.000 

45.000 

46.000 

47.000 

48.000 
49.000 10 
50.000 
51.000 


BKmon 4TH ORDBR OBSERVER S^TEM 
C0r«0N/ijt!TT/A31 , A32, A41 ,A42, 01 ,02,03,04,132 
EXTERNAL PCT.OUTP 

DIMH1SI0N y(5) , DERY(5) , PRMT(5) , AUX(1 6, 8) 

OUTPUT ' INPUT Y10,Y20' 

DATA (Y(I), I«1,4)/O.,0.,0.,0./ 

OUTPUT 'INPUT PRMT(t),I=1,4’ 

INPOT, (PRMT(l)iN,4) 

OUTPUT 'INPUT BIO' 

INPOT BIO 

A31 =>-1.622 

A32=1.1926 

A41=8.1096 

A42=-8.1096 

B32=.0437 

04=-4*BIO 

03^(G4*A31 -4*BIO**3)/A41 
02=(-BI0**4-A32*A41 •♦A42*A31 )/(-A32^»A41 tA42*A31 ) 
01=(-6*BI0^H»2-A31-A42-fA42'K>2)/(-A41 ) 

OUTPUT 01,02,03,04 
NDIM=4 

OUTPUT 'INPUT 1 TO RUN' . ‘ , 

INPUT J 

IP (J .HE. 1) 00 TO 9 
DO 1 1=1,4 
DERY(I)=0.25 

CALL HPCO(PRMT, Y, DERY, NDIM, IHLP, PCT, OUTP, AUX) 

STOP 

END 

SUBROUTINE PCT (X. Y, DERY, INO) 

DIMENSION Y(1) ,DERY(1) 

C0MM0N/PCTT/A31 , A32, A41 , A42, 01 , 02, 03, 04, B32 

C0MM0N/DAB/T1,THTD 

IP(INO .B3. 0) GO TO 2 

FORMAT (20) 

FORMAT (30) 

READ(104.4)T1,THTD 
RBAD(103,5)T2,A1.A2 
DERY(1 )=01*(-Y(4)4THTD)+Y( 

DERY(2) =02*(-Y(4) +THTD) +Y(4) 
DBRY(3)=03*(-Y(4)-tTHTD -»A31*Y(1 )-fA32*Y(2)-B32*A2 
DBRY(4)==G4*(-Y(4)4-THTD)+A41*Y(1)4A42^Y(2) 

RETURN 

END 

SUBROUTINE 0UTP(X, Y, DERY, IHLP, NDIM;PRMT) 

COMMON/DAB/T1,THTD 

DIMENSION Y(1 ),DERY(1 ),PRMT(1) 

WRITE( 1 06, 1 0)X, ?f(2) , Y(4) 'EHTD 
P0RMAT(4E13.5) 

RETURN 

END 


Balloon Fifth Order Observer 
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1.000 0 
2.000 
5.000 
4.000 
5.000 9 
6.000 

7.000 

8.000 
9.000 

10.000 

11.000 

12.000 

13.000 

14.000 

15.000 

16.000 

17.000 

18.000 
19.000 
20.Q0Q 
21.000 
22.000 
23.000 
24.000 
25.000 
26.000 
27.000 1 
28.000 

29.000 

30.000 

31.000 

32.000 

33.000 

34.000 

35.000 

36.000 4 

37.000 5 

38.000 

39.000 

40.000 2 
41.000 

42.000 

43.000 

44.000 

45.000 

46.000 

47.000 

48.000 

49.000 

50.000 
51.000 10 
52.000 
53.000 


mydm to test 5th oi^br observer system with input bias 
COr«ON/K3TT/A31 , A32,A41 >A42,GM ,02,G3f04»Q5,B32 
EXTERNAL PCT.OUTP 

DIMENSION Y(5),DERy(5).PRMT(5),AUX(l6.8) , 

OUTPUT'INPUT Y10,Y20,Y30' 

DATA (Y(l), I=1,5)/0.,0.,0.,0.,0./ 

OUTPUT ’INPUT PRMT(l),Ia1,4' 

INPOT, (PRMT(I),I*1,4) 

OUTPUT ’INPUT EIG’ 

INPUT BIG 

A31P-1.622 

A324.1926 

A41=e.1096 

A42=-8.1096 

B32~.0437 

G2=(5*EIG^4-A42*A31 +A41 *A32)/(-A42»A31 4A41 *A32) 
01=(10*EIG**2-A42^2+A31 m2)/K^\ 

G5=>-BIG*^*5/ ( A42^^A31 “A41 *A32) 
a4=-5*Ela-G5 

G3“(-1 0*1IG»*34A31 *G4+(A31 4A42) ^K>5) /A41 

OUTPUT G1,G2,G3,G4,G5 

NDIM=5 

OUTPOT'INPUT 1T0RUN’ 

INPOT J 

IP (J .NB. 1) 00 TO 9 
DO 1 I=;1,5 
DERY(I)=0.2 

CALL HPCG(PRMT,Y,DBRY,iroiM,IHiP,PCT,OUTP,AUX) 

STOP 

END 

SUBROUTINB PCT (X, Y, DERY, INO) 

DIMENSION Y(iy,DBRY(l) 

C0MM0N/PCTT/A31 , A32, A41 , A42, G1 , G2, G3, G4, 05, B32 

COMMON/DAB/TI.THTD 

IP(IN0 »EQ. 0) GO TO 2 

FORMAT ^20) 

FORMAT (30) 

READ(104,4)T1,THTD 

READ(103,5)T2,A1,A2 

DERY i)=Y(3)4G1*(-Y(4)-Y(5)+THTD) 

DERY(2)=Y(4)-HJ2*(-Y(4)-Y(5)+THTD) v 

DERY 3)=A31 *Y(1 )4A32*Y(2)4G3*(-Y(4)-Y(5)+THTD)-B32»A2 

DERYC4^4in(1 )4A42*Y(2)4G4*(-Y(4)-Yte 

DERY(5 ) =05^^ (-Y(4) -Y (5 ) +THTD) 

RETURN 

END 

SUBROUTINE 0OTP(X, Y, DERY, IHLP, NDIM, PRMT) 

COMMON/DAB/TI.THTD 

DIMMSION Y(1),DERY(1),PRMT(1) 

WRITECI 06, 1 0)X, Y(2) , Y(4) , Y(5 ) , THTD 
FORMAT (5E13. 5) 

RETURN 

END 



ORIGINAl PAGE IS 
OF POOR QUALITY 


k, Hainmlrw~Prodlcbor Corrector 


1.000 c 

2 . 0 C 0 C 

3.000 

4.000 C 
f.OOO C 

6.000 
7 . COO, 
e.ooc 
9.000 

10.000 

11.000 

12.000 

13.000 

14.000 

1 5 . 000 
16.000 

17.000 c 

18.000 C 

19.000 

20.000 
21.000 

22.000 C 

23.000 c 
24 » 000 

25.000 

26.000 C 

27.000 C 

28.000 

29.000 

30.000 

31.000 

32.000 

33.000 

34.000 0 

33.000 C 

36.000 

37.000 

38.000 C 

39.000 

40.000 

41 .000 

42.000 C 

43.000 C 

44 . 000 

45.000 

46.000 

47.000 

48.000 

49.000 

50.000 

51.000 
52.000 C 


OmiOUTINE HPCG(PPMT,y.rEPy,NriW,JHLF.FCT.OUTP,AUX) 


DIPFN 8 JOK PFKT ( 1 ) , y f 1 ) , DEPy ( 1 ) , A UX ( 1 6 , 1 ) 

N=r1 

IFIF*r.O 

X-PPMTd) 

H=PPM!T( 3 ) ‘ . 

FPWT( 5)“0 ( 

ro 1 T =1 NDIM I 

AUXf 16 .I)= 0 . \ 

AUX( 15 .l)~DFPy(l) 

1 AUX( 1 .l)sY{l) i 

IF(H*(PRMT(2)-X))3,2,4 

FEPOR RETURNi? 

? IHIFrlP 
> GOTO 4 

3 IH1F=13 

COMPUTATION OP DERY PGR STARTING VALUES 

4 IN 0:^1 

CALL PCT(X,Y,rERy,INO) 

RECORDING OF STARTING VALUES 
CALL OUTP (X , Y, DERY, IHLP, NDIM , PRMT ) 
TF(PRM.T( 5 )) 6 , 5,6 

5 JF(IHLP)7,7,6 

6 RETURN 

7 DO 8 1 ^ 1 , NDIM 

8 AUX(8,I)=riRy(l) 

COMPUTATION OF AUX( 2 .I) 

ISV =1 
GOTO 100 

9 X*X+H 

DO 10 1 - 1 , NDIM 

10 AUX(2.l)*y(l) 

JNCRFMENT H IS TESTED EY MEANS OF BISECTION 

11 IHLF=IHLFh 1 
X =X“H 

DO 12 1 = 1 , NDIM 

12 AUX( 4 , 1 )=AUX( 2 ,I) 

H=H 
N =1 
IS ¥=2 
GOTO 100 
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5^.000 

13 

X=X4H 


5^ 000 


TN0=0 - V 


55.000 


cAii pcT(x,y,rFFy.iN0‘) 


56.000 




57,000 


DO 14 I=l.NriM f 


58 000 


A<JX(2.l)=y(I) , \ 


59.000 

14 

AUX(9-I)«DIFY('I) 


60.000 


1SW=3 


61.000 


GOTO 100 


62.000 C 




63,000 C 


COMPUTATION CP TE8T VALUE CELT 


64.000 

15 

r>ELT=0. ; 


65.000 


DO 16 Isrl,NDlM 


66.000 

16 

DE1T=D11T4AUX (1 5, I)*ABS(Y(I )-AUX(4, I ) ) 


67.000 


DELTr.06666667i^DEIT 


68.000 


GO TO 19 


69.000 

17 

IF(IHLP-1 0)1 1,18,18 


70.000 C 




71.000 C 


NO SATISFACTORY ACCURACY AFTER 10 BISECTIONS. ERROR 

MESSAGE. 

72.000 

18 

IH1P=1 1 


73.000 


X=X4H 


74.000 


GOTO 4 


75.000 C 




76.000 C 


THERE IS SATISFACTORY ACCURACY AFTER LESS TEAK 11 El 

SECTIONS 

77.000 

19 

X=X4H 


78.000 


IN0=0 . 


79.000 


CALL PCT(X,Y, DEFY, INC) ' 


80.000 


DO 20 1=1, NDIM 


81 . OOC 


AUX(3. I)=Y(I) 


82.000 

20 

AUX(10.T)=DEFY(T) 


83 . 000 


N=3 


84.000 


ISW=4 


85.000 


GOTO 100 


86.000 C 




87 . 000 

21 

N=1 ■ o 


88.000 


X=X4H 


89.000 


IN0=0 


90.000 


CALL PCT(X,Y,DERY, INO) 


91 .000 


X=PFMT(1) 


92.000 


DO 22 1=1, NDIM 


93 . 000 


AUX(11,I)=DERY(1) 


94. 000 

220Y{1)=AUX(1 , I)4H*(.375«AUX(8, I)+.7916667*AUX(9,I) 


95.000 

1 

-.2083333*AUX(10,I)4.04166667*DERY(I)) 


96.000 

23 X=rX4H 


97.000 


N=N4l 


98.000 


IN0=1 


99. 000 


CALL EOT (X,Y, DEFY, INO) 


100.000 


CALL OUTP(X, y, DEFY, IHLP,NDIM, PRMT) 


101 .000 


IF(PPMT(5))6,24,6 


102.000 

24 

IF(N-4)25, 200, 200 


105.000 

25 

DO 26 1=1. NDIM 


104.000 


AUX(N,l)=y(l) 


105. OCO 

26 

AUX(N47, I)=DEFY(I) 


106.000 


lF(N-3)27.29*200 



109 


107.OCO C 
108.000 

109.000 

110.000 
111.000 
112.000 

115.000 C 

114.000 

115.000 

116.000 

117.000 

118.000 

119.000 C 

120.000 C 

1 21 . 000 C 

122.000 C 

125.000 

124.000 

125.000 
126.000 

127.000 C 

126.000 C 

129.000 

150.000 

151.000 

152.000 

155.000 

154.000 

155.000 

156.000 C 

157.000 

158.000 

159.000 

140.000 

141.000 

142.000 

145.000 

144.000 C 

145.000 

146.000 

147.000 

148.000 

149.000 

150.000 

151 .000 

152.000 C 

155.000 C 

154.000 C 

155.000 C 

156.000 C 

157.000 

1 58.000 

159.000 C 

160.000 C 
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27D0 28 I*1.NDIM OF POOR QUALITY 

DEIT UX ( 9 , 1 ) -f A UX ( 9/1 ) /p 

DEXTsDILT+DEIT 

28 y(I)xA 8 X( 1 ,l)+. 5555535^H»(AUX(e,l)+DEIT4 AUX( 10 ,I)) 
GOTO 25 


29 ro 50 Ie1,lIDIM 

PEIT^AUX ( 9 , 1 ) -fA UX ( 1 0 , 1 ) 

DFIT=DFXT+rEIT4rEIT 

50 y(I)-AUX(1 , 1)4. 575^fH*(ArX(e. I)4I'F1T4A0X(1 1 , X) ) 
GOTO 25 


;\ 


THE FCLIOVIXG PAPT OF SUPPOUTINF HPCG COHFUTFS BY MEANS OF t 
PliHGE-KlTTTA MFTHOP STARTING VAIUES FOR THE NOT SELF- STAPTINQ 
PPmCTOR* CORRECTOR MFTHOr. I 

100 DO 101 I=1,NriM I 

Z=H*AUX(N4?,I) 

AUX{5.I)=Z 

101 Y(I)=rAUX(N. I)4.4*2 

Z IS AN AUXILlAPy STORAGE LOCATION 

Z"X4.4*H 

INO^O 

CALI FCT(Z,Y, DERY, INC) ^ 

DO 102 1 = 1 ,NDIM ^ 

Z=H*DEFy(l) . 

AUX( 6,i‘)=Z 

102 Y(I)=AUX(N,l) 4 . 2969776 #AUX{ 5 ,l) 4.15e75964Z 

Z=X4.4557572*H 

IN0=0 

CALL PCT(Z,Y, DEFY, INC) 

DO 105 I=1,NDIM 
Z=H*DERY(l) 

AUX(7,I)=Z ‘ 

105 yfI)=AUX(N,I)4.2181004*AUX(5,l)-5.050965*Al‘X(6,I)45.852865<^^ 

Z=X4H 

IN0=0 

CAIL FCT(Z,y,DERy,INO) 

DO 104 I=1,NPIM 

1O4Oy(I)=AUX(N,I)4.1747605*ArX(5. l)-.5514e(&7«AUX(6. 1 ) 

141 .20555 6 #AGX{ 7 .l) 4.1711846»H^»DEPy(l) r 

G0T0(9, 15- 15,21 ).ISW 

POSSIBIE BREAK-POINT FOR IINKAGI / 

STARTING VAIUES ARE COMPUTED. i/ 

NOW START HAMMINGS MODIFIED PREDICTOR'^- CORRECTOR METHOD. 

200 ISTFP=3 

201 IF(N-8)2C4,202,204 ’ 

N=8 CAUSES THE ROWS OF AUX Tp CHANGE TKElR STORAGE L0CAT10Ns| 
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161 .000 
162.000 
163.000 
164.000 
163.000 
166.000 
167.000 
166.000 
169.000 
170.000 
171.000 
172.000 

173.000 

174.000 
175.000 
176.000 

177.000 

178.000 
179.000 
160.000 
161.000 
182.000 

183.000 

1 84 . 000 

1 85 . 000 
186.000 
187.000 
188. COG 
188.000 
190.000 
191 .000 
192.000 
1 93 . 000 
1 94 . 000 

1 95 . 000 

196.000 
1 97 . 000 
1 98.000 
1 99.000 
200.000 
201.000 
202.000 
203.000 

204 . 000 

205.000 
206.000 
207.000 
208.000 
209.000 
210.000 
211.000 
212.000 
213.000 
214.000 


ORKSiNAt 

202 DO 203 Ns2,7 OF POOR Q«Ati ^ 

DO 203 I“1,NDIf« 

AUX(N-M)s^Al»X{N,I) 

203 AUXCjJ-i6, 1)~ADX(N47.t) 

N*t7 

C 

C N IE8S THAN 8 CAD6E8 N+1 TO GET N 

204 N:^N4l 
C 

C COMPUTATION OP NEXT VECTOR Y 
DO 205 I=1,NDIM 
AUX(N-1,1 *Y(I) 

205 AUX(N46,l)=DERy(l) 

X=X-fH 

206 ISTEP=ISTEP4l 
DO 207 I=1,HDIM 

0DELT-AUX{N-4, 1)4l .333333*H4(AUX{N46, I)4AUX(N46,I)-AUX(N45,I) 
1AUX(N44,T)4AUX(N44,1)) 

Y ( I ) =DEXT- .92561 98*A UX ( 1 6, 1 ) 

207 AUX(16, I)-DEIT 

C PREDICTOR IS NOW GENERATED IN ROW 16 OF AUX> MODIFIED PREDIC 

C IB GENERATED IN Y. DELT MEANS AN AUXIEIARY STORAGE. 

C 

IN0=0 

CAIL PCT(X,Y,DEPY,INO) 

C DERIVATIVE OP MODIFIED PREDICTOR IS GENERATED IN DERY 
C 

DO 208 I=1,NDIM 

ODEIT=.125*(9.*AUX(N-1. I)-AUX(N-3,I)43. *H*(DEPYfl)4AUX(N46.l)- 
1 AUX(Nh6,I)-AUX(N45.I))) 

AUX(16,l)=rAUX(16,l)-DFIT 
206 Y ( I ) =DFLT 4 . 074 3 80 1 7 ^^A UX ( 1 6 . 1 ) 

C 

C TEST WHETHER H MUST BE HALVED OR DOUBLED 
DELT-0. 

DO 209 I=1,NDIM 

209 DELT=DELT4A UX ( 1 5 , 1 ) *ABS ( A UX (16.1)) 

GO TO 210 

C • 

C H MUST NOT BE HALVED. THAT MEANS Y{l) ARE GOOD. 

210 IN0^1 ^ 

CALL PCT(X,Y,DERY, INO) 

CALL OUTP (X, Y, DERY, IHIP, NDIM, PRMT) 

IP(PRMT(5))212,211,212 

211 IF(IHLP-1 1 )213, 212,212 

212 RETURN 

213 IF(H*(X-PRMT(2)))214,212,212 

. 214 IP(AES(X-PRMT(2))-.1*AES(H))212,215,215 
215 GO TO 201 

C : . . . ■ 

C 

C H COULD BE DOUBLED IP ALL NECESSARY PRECEEDING VALUES APE 
C •AVAILABIE 

216 IP(IHLP)201, 201, 21 7 


lit 


215-000 

216.000 

217.000 

218.000 

219.000 

220.000 
221.000 
222.000 
225.000 
224 . OOC 
225- 000 
226.000 

227 . 000 

228.000 
229-000 
250.000 
251 *000 
252.000 

255.000 

254.000 C 

255.000 C 

256.000 G 

257.000 

258.000 

259.000 

240.000 

24 1 . 000 
24 2.000 
245-000 

244.000 
245-000 

246.000 

247.000 

248.000 
249-000 

250.000 

251 .000 

252.000 

255.000 

254 .000 

255-000 

256.000 

257.000 

258.000 

259.000 

260.000 
261 ,000 
262.000 
265-000 
264 - 000 
265.000 
266.000 
267.000 
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217 IF(F-7)201 .218,216 ’ ' ^ OT POOR QUALITY 

218 IF(J8TFr.4)201,219-219 

219 IM 0 r . I 8 TFF /2 

IFf T8TFP- IFCr- IMCr)20l . 220. 201 

220 H*H . f \ 

IHiFsrJHlF-1 .i 

T8TFF=0 / 

PC 221 T:r1,FFTM 
AUX(H-1,I)“AltX(N 2,1) 

AUX(N.2. 1 =AliX(N 4,1 . 

AGXfl<-5.l5'-AUX(II-6,I) f 

AI>X(N+6,T)=AUX(N45.I) 

AUX(F45.I =rAUX(K+5.I 

A‘UX(N44,l)=rAUX(N4l,I) ^ ,( 

DELfsA UX (N46 . 1 ) 4A UX (N45 , 1 ) 

DEIT -EEIT 4DELT 4DEIT 

2210AlTX(l6.I)-8.962965l'(y{I)-AUX(N..5,l))-5.361111*H*(DERY(J)4DEIT 

14AGX(N44.I)) 

GOTO 201 


H MU8T EF EAIVEE 
222 JHLP=IHIP41 

TP(IFLF-10)225, 225*210 
225 HjtF 
T8TFP=0 

DO 224 J==1.NDTF • * 

Oy(X)^. 00590625*^80. 1^AUX(N-1, 1)41 55 . -^AUXOlJ-a, 1)440- *AUX(N-5, 1)4 

I)-AUX(N44, 1))*H , 

0AITX(F-4,I)='.00590625*(12.*AUX(N-1.I)4155.*AUX(F-2,I)4 
1108. *A UX f N- 5 , 1 ) -1 A UX (N-4 , 1 ) ) . 0254375* ( A UX (N46 . 1 ) 41 8 . *A UX (N45 .1)- 
29.*AUX(N44. T))*H 
AUX(N.5.I)*AITX(F«2.I) 

224 AUXfN44.T)=AlOr(N45,I) 

X:^X-H 

DFIT=X-(H4H) 

IN0=0 

CALI FCT(DELT.y,DEFy INO) 

DO 225 I=1.NDIF I 

AUX(F-2,I)=Y(]) ’I 

AUX(N45. l)=UERY(l) I 

225 Y(I;=AUX(F- 4 ,I) 

DEIT=DEIT.. (H4H) 

iNo=o V ' r 

CALL FCT(DELT,y,DERY,INO) I 

DO 226 I=1,NDIM I 

DELT =A UX (F 45 , 1 ) 4 A UX (N4 4 , 1 ) 

DELT=^DELT4DELT4DFLT I 

0AUX(16,I)=8.962963*(AUX(N-1,I)-Y(I))-5.361111 W(AUX(N46,I)4DEL'^^ 
l4DFRy(I)) I 

226 AUXfF45,I)=DPPyM) 

GOTO 206 

FFD ’ 


. ORIGINAL PAGE 18 

^ ' ^ OF POOR QUALITY 
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